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SUMMARY 


This  report  details  the  way  a  multistate  ref lectometer  is 
used  to  measure  the  reflection  coefficient  of  waveguide  devices. 
It  contains  theoretical  derivations  of  the  equations  and  a 
computer  program,  with  explanatory  notes,  for  solving  them. 

The  calibration  procedure  is  also  given.  These  techniques 
are  suitable  for  all  frequency  bands  and  transmission  media 
where  suitable  standards  and  detectors  exist. 
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1  INTRODUCTION 

(1  2) 

The  theory  of  multistate  ref lectometer s  has  been  well  documented 
This  paper  demonstrates  the  practical  Implementation  of  that  theory.  The 
circuit  arrangement  is  a  recent  development  suggested  by  G  F  Engen,  and  has 
the  advantage  over  earlier  designs  of  a  wide  range  of  calibration  techniques. 
The  theory  of  this  circuit  Is  given  and  some  of  the  Implications  discussed. 
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The  rest  of  the  paper  shows  the  equipment  and  computer  program  used  to 
construct  a  multistate  ref lectometer  for  measuring  the  voltage  reflection 
coefficients  of  waveguide  devices.  The  program  for  operating  the  ref lecto¬ 
meter  and  performing  numerical  operations  on  the  data  is  entitled  "MULTI STATE" 
version  2  and  has  been  written  in  Basic  4.0  for  use  with  the  4000  series 
Commodore  PET. 

2  MULTI STATE  REFLECTOMETER  THEORY 

Figure  1  illustrates  the  multistate  ref lectometer . 


Assume  that  each  position  k  of  the  reflector  in  Figure  1  produces  a 
different  repeatable  state  of  linear  transmission  and  reflection  within  the 
waveguide  junction  having  ports  1  to  4.  Postulating  that  at  any  port  there 
is  no  evanescent  mode  and  only  one  transmitted  mode  (but  not  necessarily  the 
same  transmitted  mode  at  all  ports),  enables  each  linear  state  to  be  described 
by  a  scattering  matrix  Ts]^.  Then,  if  a^  and  b^  respectively  represent  waves 
incident  on  and  emergent  from  port  i  (where  i  ■  1,2, 3, 4);  [b]  ■  taj^t*]. 
Denoting  the  reflection  coefficients  of  the  detectors  providing  power 
indications  P^  and  P^  by  ■  *3^3  and  ■  a ^/b^,  and  substituting  gives: 
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Since  there  is  non-identical  transmission  between  the  parts  1  to  4,  the 
square  matrix  may  be  inverted  for  each  k.  Denoting  the  inverse  by  [M]  and 
its  elements  by  m^: 
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Thus  for  the  multistate  4-port  junction  1,2, 3, 4  for  each  state  k  (suppressing 
the  suffix  k  for  the  moment) 


'i  “  milbl  *  mi2b2 


bi  “  milbl  *  mi2b2 


(for  i  =  1,2) 

(where,  now,  i  -  3,4) 


(2.3) 


whence 


i  -  3,4) 


..  mil  /  m22mil\  . 
bi  “  5^  *2  +  y  i2  -  S^-J  b2  (for 

Substituting  T  ■  a2^b2  8ive® 

(for  i  «  3,4) 


(2.4) 


bi  "  ^®ir  *  ®i^b2 


(2.5) 


where  j  ^"*21  *nd  ®i  “  *mi2  “  m22mil^m21^ 

Reintroducing  the  suffix  k  gives 
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(2.6) 


(Vk  “  (V  4  Vk  b2 


Similar  analysis  of  the  invariant  3-port  junction  2,3,5  produces: 


(<*3r  +  0j)  b2 


(2.7) 


where,  now,  et^  and  Bj  are  dependent  on  only  the  scattering  coefficients  of  the 
3-port  junction  and  the  detector  at  port  3.  Combining  equation  (2.7)  and 
(2.6)  gives 


(Vk 


V  +  ek 


(2.8) 


Equation  (2.8)  also  describes  a  ref lectometer  with  (k+1)  detectors  connected 
to  a  time-invariant  junction  having  (k+3)  ports.  For  k*3  the  six-to-four 
port  reduction  algorithm  can  be  used,  making  calibration  with  three  precisely 


known  loads  applicable 


(3,4) 


Those  methods  of  calibration  were  not  chosen  in  this  case  as  speed  of 
calibration  was  not  a  criterion  and  the  extra  redundancy  in  the  method  chosen 
gives  a  better  assessment  of  the  uncertainties  involved  in  calibration  and 
measurement. 
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Substituting  peJ  -  T  and  converting  from  complex  to  polar  notation 
equation  (2.8)  becomes: 


/PA  A  *  (PPJZ  ♦  2P 

(4  '  s  {—SsTZT, 


2pR^  cos(Pj{'f0)' 


(2.9) 


S  cos(F+0) 


It  is  useful  from  the  viewpoint  of  calculating  uncertainties  to  assume  that  S 
and  F  are  not  invariant  but  are  dependent  upon  the  state  thus: 


fa\  m  L  tl  ♦  (pV*  4  2pRk  co,(Ve)\ 

V  3/k  k  ll  ♦  (pSfc)*  ♦  2pSk  cos(Fk+0y 


(2.10) 
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This  is  s  non-linear  equation  with  5  unknown  L,R,P,S,F  for  each  state  k. 

A  solution  could  be  found  immediately  using  an  iterative  least-squares  fitting 
method.  But  a  more  accurate  method,  if  low  values  of  reflection  coefficient 
are  to  be  measured,  is  to  determine  the  k  values  of  L  and  then  use  the  fitting 
technique  to  solve  for  each  R,P,S  and  F. 

2.1  FINDING  VALUES  OF  L 

With  a  perfect  load  on  the  output  port  p*0  and  suppressing  the 
suffix  k  again,  equation  (2.10)  becomes  P4/P3  *  L.  Finding  a  perfect 
load  is  difficult  so  one  is  emulated  using  a  fairly  good  load  and 
precision  spacers. 

If  p  is  small  and  the  directivity  of  the  second  coupler  is  good 
(giving  S  <  .05)  then  the  term  (pS)  becomes  negligible  giving 

Pa  2 

—  *  L  (1  ♦  (pR)  2pR  cos(P+0))  (2.11) 

p3 

Then  introducing  known  phase  shifts  a  and  0  by  inserting,  in  turn,  two 
precision  spacers  between  the  output  port  and  the  load  gives  3  equations 
to  solve: 


L  (1  +  (pR)2  +  2pR  cos(P+0)) 


L  (1  +  (pR)2  +  2pR  cos(P+0+a)) 


L  (1  +  (pR)2  ♦  2pR  cos(P«-e+0)) 


(2.12a) 


(2.12b) 


(2.12c) 


Let  Sj  -  (p^/PjJj.  «2  *  (p4/p3>2  *nd  *3  “  ^p*/p3*3  *nd  comblne  p 
and  6  into  a  single  phase  shift  6',  then  let 
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cos6'  -  cos(8'+ot) 
cos6 '  -  cos(6'+B) 


(2.13) 


C1 


sl~  s2 
S1 '  83 


c^(l-cosB+tan9' sinB)  *  l  -  coso +  tan9' sina 


(2.14) 


tan©1 


c^(l-cosfJ)  +  cosa-  1 
sina  -  c^(sinS) 


(2.15) 


From  this  the  value  of  9’  is  found.  If  we  then  let  C2  *  s^/s^  then 


1  +  (pR)  +  2pR  cos9' 

.  ,  - - 

1  +  (pR)  +  2pR  cos(0'+0) 


(2.16) 


Let  Cj  *  C2  cos(9'+0)  -  cos 9 1  then 


(cj-lMpR)  +  2c^pR  +  c2“^  *  0 


(2.17) 


This  is  a  simple  quadratic  in  pR  and  by  putting  c^  *  2cj/(c2-1)  the 
required  solution  is 


(2.18) 


knowing  9'  and  pR  then  L  is  calculated  from  equation  (2.12a). 


2.2  FINDING  VALUES  OF  R,P,S  AND  F 


Having  determined  L  a  movable  short  is  used  to  evaluate  the  remaining 
4  parameters.  Suitable  redundancy  for  the  iterative  least  squares  fitting 
technique  is  provided  by  using  7  data  points  for  each  state.  To  minimise 
the  number  of  iterations  required  to  reach  the  optimum  solution  a  good 
set  of  starting  estimates  are  required.  These  are  called  the  starting 
parameters  and  a  new  set  are  required  for  each  state. 
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2.2.1  Finding  Starting  Parameters 


This  part  of  the  calibration  procedure  is  carried  out  using  a 
short  circuit  with  various  phase  shifts  supplied  using  spacers  £>r 
a  movable  short  circuit)  so  p  is  replaced  by  -1.  Let  Y  *=  p^/p^  *0T 
a  particular  phase  shift  Xt  then  equation  (2.10)  becomes 


2R  cos (P+X) 


S  -  2s  cos(F+X) 


) 


f  2 .19) 


A  good  set  of  starting  parameters  can  be  found  by  treating  equation 

2 

(2.19)  as  a  linear  equation,  ie  taking  R  and  R  as  independent 
parameters . 


Let  a  *  1+R  ,  b  »  1+S  ,  c  = 

and  d  *  -2R  sinP  ,  e  =  -2S  cos  F 


-2R  cosP  , 

and  f  =  -2S  sinF 


Then  equation  (2.19)  becomes 


H 


a  +  c  cosX  -  d  sinX\ 
b  +  e  cosX  -  f  sinX/ 


(2.20) 


j-  *  si  +  s2  cosX  -  s3sinX  -  a1*  y-  cosX  +  s5  sinX  (2.21) 

L  L  L 

where  si  a  a/b,  s2  -  c/b,  s3  «  d/b,  s4  =  e/b  and  s5  =  f/b. 

A  fifth  order  simultaneous  equation  solver  (using  elimination 
and  having  no  safeguards  against  the  matrix  being  singular  or 
making  any  attempt  to  optimise  the  sequence  of  elimination)  is 
used  to  evaluate  si,  s2  etc  from  equation  (2.21).  From  the  original 
substitutions 


P  -  tan”^(s3/s2)  and  F  »  tan"^(s5/s4) 


( R/S } 


a22+s32 

s42+s52 


and 


1  +  R 
1  +  S2 


(2.22) 
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Then 


2.2.2  Least  Squares  Fitting  Procedure 


The  least  squares  fitting  is  done  by  taking  a  Taylor  expansion 
about  each  of  the  parameters  thus  linearising  the  equation  to  enable 
normal  equations  to  be  written  that  can  be  solved  to  obtain 
corrections  to  the  original  estimates  of  the  parameters.  To 
implement  this  method  the  partial  derivatives  with  respect  to  each 
of  the  parameters  are  required. 


From  equation  (2.19) 


3Y 

2L(R-cos(P+X)) 

m  * 

1+  S2  -  2S  cos(F+X) 

3Y 

2LR  sin(P+X) 

5p  “ 

0 

1 +  S  -  2S  cos(F+X) 

8Y 

-2Y(S  -  cos(F+X)) 

3s  “ 

1  +  S  -  2S  cos(F+X) 

8Y 

-2YS  sin(F+X) 

5F  ” 

1 +  S2-  2S  cos(F+X) 

(2.24a) 


(2.24b) 


(2.24c) 


(2.24d) 


t  \ 

I 

i 


In  terms  of  the  estimates  of  the  parameters  (making  allowance  for 
them  being  incorrect)  we  can  write 


Y 


1  4 
.1  + 


R  -  2R  cos(P+X) 
- - 

S  -  2S  cos(F+X) 


:]*d) 


* ♦  (f ) » * 


(2.25) 


r 


•4 
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m  -  Y  -  L 


R  -  2R  cos (P+X) 
S2  -  2S  cos (F+X) 


(S) « *  (S) » *  (S)  -  *  (¥) 


(2.26) 


For  a  given  estimate  of  the  parameters,  m,(8Y/8R)  etc  are  all  known 
and  8R  etc  are  unknown. 

The  least  squares  method  applied  to  equation  (2.26)  proceeds 
by  minimising 


i*NP 

IX 


where  NP  is  the  number  of  points  and 


[. . -l2  .[ . 

L  J Point2  L 


Pointl 


(2.27) 


This  means  we  are  searching  for  a  stationary  (minimum)  value  of  I  e^ 
given  by  a  particular  set  of  6R,  6P,  6S,  and  6F.  By  taking 


die. 


■0T6RT  "  °’  etc  • 


we  obtain  the  normal  equations: 
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[S-lH(i)2]  “♦[S-S]**[S-5]*r 

[i-n]«*[lf-l]«*[f-s]-[(f)2]- 

where  the  []  brackets  indicate  summation  over  all 


data  points. 


The  four  simultaneous  equations  are  solved  for  6R,  6S ,  6P  and 
6F  using  the  determinant  method,  which  are  applied  as  corrections 
to  the  last  estimates  of  R,  S,  P  and  F  used.  The  sum  squared  error 
is  then  checked  and  if  the  change  from  the  last  cycle  is  <  5  x  10 
then  the  current  values  of  R,  S,  P  and  F  are  accepted  as  the  best 
estimate . 


2.3  STANDARD  DEVIATION  OF  THE  INDIVIDUAL  PARAMETERS 


If  RQ  is  the  total  RMS  error  of  the  points  then  the  standard 
deviation  of  the  parameters  is  given  by  the  product  of  RQ  and  the  ratio 
of  two  determinants  D(P)  and  D2^  .  D2  is  the  4x4  determinant  used  in 

the  equation  solver  and  D(P)  is  a  3x3  determinant  created  from  D2  by 
removing  the  relevant  row  and  column.  Thus  D(R)  is  D2  with  the  first 
row  and  column  removed  etc. 

This  cycle  of  iterative  least  squares  fitting  and  then  finding  the 
standard  deviation  of  the  individual  parameters  is  performed  for  each 
state. 
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2.4  MEASUREMENT  OF  AN  UNKNOWN 


Reintroducing  the  suffix  k  to  denote  the  state,  then 


Yk  1+  (pRk)  +  2pRk  cos(Pk+6) 
**  Lk  1+  (pSk)2  + 2pSk  cos(Fk+9) 


(2.29) 


'  where  p  and  6  are  now  the  unknowns. 


2.4.1  Starting  Parameters 


We  linearise  the  equation  by  considering  (pR  )  as  independent 

2  K 
of  pRj^  and  (pSk)  as  independent  of  pSk>  Thus  expanding  equation 

(2.29)  for  3  states  we  obtain 


2  2 

m^-1  *=  sl(Rj-mjS^)-fs2(R^cos(P^)-m^S^cos(F^))  + 


s3(m1S1sin(F1)-R1sin(P1)) 


2  2 

n^-l  *  8l(R2-m2S2)+s2(R2Cos(P2)-m2S2C08(F2))  ♦ 


s3(m2S28in(F2  l-RjSin^ )) 


2  2 

m^-l  ■  sl(R^-m2S2)+82(R2Cos(Pj)-in^S2Cos(F2))  + 


(2.30) 


sS^jS^sin^J-R^sin^) ) 


where  il  ■  s2  -  2pcos0  and  s3  «  2psin6. 

Hence  starting  values  required  are  VsT  and  tan”*  s3/s2. 


The  3  simultaneous  equations  are  solved  by  the  determinant 
method . 
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Figure  2  shows  the  general  layout  and  equipment  used  for  this  multistate 
ref lectometer .  The  operation  of  most  of  the  equipment  is  implicit  in  its 
function  but  a  more  detailed  look  at  some  of  the  non-standard  components  is 
helpful. 
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3.1  THE  BENTHAM  SMD-3D  DRIVER  AND  STEPPING  MOTORS 


The  motors  used  have  12  Volts  4  phase  windings  and  give  a  rotation 
of  1.8°/step.  The  driver  has  a  front  panel  switch  giving  capability  of 
selecting  speeds  between  1  step/second  and  400  steps/second.  This 
facility  is  not  fully  available  via  the  IEEE-488  bus.  Only  one  of  two 
internally  settable  speeds  can  be  selected.  The  instructions  to  the 
driver  must  be  of  the  form  "XNNNNY",  where  X  is  either  "D"  or  "I" 
denoting  direction  and  Y  can  be  any  one  of  four  characters  indicating 
which  motor  and  which  preset  speed.  The  "NNNN"  signifies  any  4  digit 
number  in  the  range  "0000"  to  "9999",  which  denotes  the  number  of  steps 
to  be  taken. 

3.1.1  Stepping  Motor  1 

This  is  coupled  to  a  waveguide  switch  via  a  general  drive 
mechanism.  This  ensures  rapid  transitions  and  good  alignment 
with  minimum  mechanical  shock. 

3.1.2  Stepping  Motor  2 

This  is  used  to  control  the  position  of  the  side-arm  short 
circuit.  A  1:1  gearbox  with  nylon  gears  is  used  to  reduce  noise 
and  allow  the  micrometer  to  move  in  and  out  without  having  to 
mount  the  motor  on  a  movable  base.  The  micrometer  thread  is 
2  turns/mm  giving  a  resolution  of  0.0025  mm/step.  The  resettability 
is  considerably  better  than  that,  when  precautions  are  taken  to 
minimise  backlash.  The  present  arrangement  gives  a  maximum  speed 
of  adjustment  of  the  short  circuit  of  1  mm/second.  Increasing  this 
speed  would  degrade  the  resolution. 

3.2  MICROWAVE  FREQUENCY  STABILISATION 

The  microwave  sources  used  are  Systron-Donner  YIG  tuned  Gunn  diode 
devices.  They,  are  voltage  tunable  across  a  waveguide  band  and  have  a 
separate  coll  for  frequency  modulation.  The  EIP  578  source  locking 
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counter  has  a  COARSE  TUNE  output  giving  0  to  4-10  V  and  a  FINE  LOCK 
output.  If  the  COARSE  TUNE  output  supplies  the  tuning  voltage  of  the 
source  and  a  voltage-to-current  converter  is  connected  between  the 
FINE  LOCK  output  and  the  FM  coil,  then  frequency  locking  can  be  readily 
achieved.  This  system  gives  excellent  frequency  stability  and  the 
facility  of  computer  selection  of  any  frequency  in  the  band. 

4  DETAILS  AND  EXPLANATION  OF  COMPUTER  PROGRAM 


Lines  100-140 

Lines  200-230 

Line  240 

Line  250 


Program  name  and  details. 

Assign  channels  for  IEEE-488  bus. 

Send  voltmeter  instructions  for  range,  triggering  and  SRQ. 

Set  up  arrays.  In  this  version  of  Basic  arrays  of  less  than 
10  elements  do  not  need  to  be  declared. 


Lines  270-320  Input  data,  frequency  and  waveguide  band  information. 

Determine  relevant  waveguide  broad  dimension.  This  version 
covers  3  waveguide  bands  from  8.2-26.5  GHz. 

Line  330  Calculate  phase-change/mm  for  that  frequency  in  that  waveguide. 

Lines  340-380  Print  headings. 


Line  390  Decide  to  perform  new  calibration  or  use  constants  stored  on 

a  disk  from  a  previous  calibration  at  the  same  frequency. 

Lines  410-550  USING  A  PREVIOUS  CALIBRATION  STORED  ON  DISK. 


Lines  430-460 


Generate  file  name  and  open  channel  to  disk  drive. 
Ensure  file  exists  and  check  for  any  disk  errors. 


Lines  470-490 


Acquire  micrometer  settings  of  the  side-arm  short-circuit  for 
each  state  and  the  calibration  constants  at  each  setting. 
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Lines  500-550 


Print  headings  and  calibration  data 


Lines  570-2300 


Lines  590-610 


Lines  630-920 


Lines  650-680 


Lines  690-820 


Lines  830-920 


Lines  940-2180 


Lines  960-980 


Lines  1000-1010 

Lines  1020-1120 


Line  1160 

Lines  1180-1540 

Line  1200 


NEW  CALIBRATION. 

Input  required  micrometer  settings  and  print  them. 

MEASUREMENT  OF  L  (section  2.1). 

Input  length  of  spacer  inserted  between  the  output  port  and 
the  attached  load  and  measure  P^/P^  for  each  state  using 
subroutine  one.  Repeat  for  two  further  spacer  lengths. 

Use  circle  fitting  technique  to  find  the  value  of  L  for  each 
state. 

Print  headings,  data  and  final  values  of  L. 

DERIVE  VALUES  OF  R,  P,  S  and  F  (section  2.2). 

Input  information  about  the  short  circuit  used  for 
calibration  and  the  number  of  data  points  required. 

Using  seven  points  seems  to  give  good  results  but  any 
number  between  5  and  20  allowed. 

Print  headings. 

Measure  P^/P^  *n  «ach  state  using  subroutine  one  for  each 
position  of  the  calibration  short  circuit.  Reverse  the  order 
of  the  states  on  alternate  cycles,  and  print  accumulated 
data. 

Derive  calibration  constants  for  each  state  in  turn. 

FIND  STARTING  PARAMETERS  (section  2.2.1). 

Select  data  related  to  this  state. 


15 


Lines  1210-1230  Set  up  matrix  of  5  simultaneous  equations.  Using 
equation  (2.21)  and  the  first  five  data  points. 

Lines  1240-1540  Use  determinant  method  to  extract  initial  values  of  R,  P, 

S  and  F. 

Lines  1560-2010  FOUR  PARAMETER  ITERATIVE  LEAST  SQUARES  FITTING  (section  2.2.2). 

Line  1580  Set  number  of  data  points  and  criterion  for  successful 

convergence. 

Line  1590  Perform  up  to  nine  Iterations. 

Lines  1620-1730  Find  partial  derivatives  at  each  data  point  and  sum  over 
all  data  points  (equation  (2.28)). 

Lines  1740-1910  Find  corrections  si,  s2,  s3,  s4  to  the  current  values  of  the 
parameters  using  the  determinant  method  of  simultaneous 
equation  solving. 

Line  1920  Apply  corrections  to  current  values  of  parameters. 

Lines  1930-1970  Calculate  the  RMS  value  of  the  differences  between  predicted 
and  measured  values. 

Line  1980  Check  if  either  the  absolute  value  of  the  RMS  deviation  or 

the  change  in  RMS  deviation  between  this  iteration  and  the 
previous  one  are  within  the  convergence  criterion.  If 
converged  drop  out  of  loop  to  Line  2070. 

Lines  1990-2010  Print  latest  values  and  go  back  to  1590  for  next  iteration. 

Lines  2050-2060  Print  warning  if  number  of  Iterations  reaches  9  without  the 

convergence  criterion  being  satisfied. 

Lines  2070-2100  Ensure  R  and  S  are  positive,  and  P  and  F  lie  between  0  and 

2s. 


Lines  2110-2150  Calculate  the  individual  standard  deviations  of  the 
parameters  (section  2.3). 

Lines  2160-2170  Store  calibration  data  for  this  state. 

Line  2180  Go  back  to  1200  until  all  states  calibrated. 


Lines  2190-2300  Print  headings  and  values  of  the  calibration  constants  for 
each  state. 


Lines  2310-2320  If  required  store  calibration  data  and  constants  on  disk, 
using  subroutine  eight. 


Lines  2340-3100 

Line  2360 

Lines  2370-2380 

Line  2390 

Lines  2410-2620 


MEASUREMENT  OF  UNKNOWN  (section  2.4). 

Set  number  of  parameters  to  2  and  the  number  of  data  points 
to  the  number  of  states. 

Input  identity  of  unknown  and  print  headings. 

Measure  P^/P^  for  e«ch  state  using  subroutine  one. 

FIND  STARTING  PARAMETERS  (section  2.4.1). 


Lines  2430-2520  Set  up  matrix  of  three  simultaneous  equations  using 
equation  (2.30)  and  the  first  three  data  points. 

Lines  2560-2600  Solve  the  simultaneous  equations  using  the  determinant 
method. 

Lines  2610-2620  Derive  starting  values  for  UM  and  UP  where  UM  and  UP  are 
the  unknown  magnitude  and  phase  respectively. 

Lines  2640-2890  TWO  PARAMETER  ITERATIVE  LEAST  SQUARES  FITTING  (section  2.4.2). 


Line  2660 


Set  number  of  data  points  and  criterion  for  successful  convergence. 


Perform  up  to  nine  Iterations. 
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Line  2670 


Lines  2680-2760  Find  partial  derivatives  at  each  data  point  and  sum  over 
all  data  points  (equation  (2.34)). 

Lines  2800-2810  Solve  2  simultaneous  equations  and  apply  corrections  to 
current  values  of  the  parameters. 

Lines  2820-2850  Calculate  the  RMS  difference  between  predicted  and  measured 
values . 

Line  2860  Check  convergence.  If  converged  drop  out  of  loop  to  Line  2920. 

Lines  2870-2890  Print  current  values  of  parameters  and  RMS  Deviation  and 
return  to  2670  for  next  iteration. 

Lines  2900-29100  Print  warning  if  number  of  iterations  reaches  9  without 
convergence  criterion  being  satisfied. 

Lines  2920-2930  Ensure  |p|  of  unknown  (UM)  is  positive  and  phase  of  unknown 
(UP)  lies  between  0  and  2«. 

Lines  2940-3080  Print  out  results. 

Lines  3100-3110  Either  go  to  Line  2340  if  another  measurement  is  required 
or  else  close  all  channels  and  program. 

Lines  3180-4040  SUBROUTINES . 

Lines  3180-3370  SUBROUTINE  ONE. 

Line  3200  Measure  the  bridge  voltages  with  no  microwave  power  present 

using  subroutine  two. 

Line  3210  Store  bridge  voltages. 

Line  3220  Turn  microwave  power  on  using  subroutine  five. 


Lines  3230-3290  Measure  the  bridge  voltage  with  microwave  power  present. 

Store  these  voltages.  Advance  side-arm  short  circuit  to 
next  position  using  subroutine  six.  Continue  this  loop 
until  bridge  voltages  measured  and  stored  for  every  state. 

Line  3300  Turn  microwave  power  off  using  subroutine  five. 

Lines  3310-3320  Measure  bridge  voltages  in  the  absence  of  microwave  power 

using  subroutine  two  and  store  them. 

Lines  3330-3360  Calculate  the  amount  of  microwave  power  indicated  by  each 
bridge  in  each  state  assuming  a  linear  drift  with  time 
between  each  pair  of  "power  off"  readings.  Find  the  ratio 
of  indicated  powers  for  each  state  and  store  them. 

Line  3370  Return  to  main  program. 

Lines  3390-3470  SUBROUTINE  TWO. 

Line  3410  Open  all  channels  on  data  scanner. 

Line  3420  Determine  the  number  of  voltage  integrations  to  be  performed. 

Lines  3430-3440  Measure  the  bridge  voltage  of  the  forward  going  power  sensor 
and  store  it. 

Lines  3450-3460  Likewise  for  the  reverse  power  sensor. 

Line  3470  Return  to  subroutine  one. 

Lines  3490-3550  SUBROUTINE  THREE. 

Line  3510  Place  voltmeter  in  averaging  mode  and  wait  until  voltages 

have  settled. 

Line  3520  Trigger  voltmeter  and  wait  for  SRQ  using  subroutine  four 

to  Indicate  completion  of  single  measurement.  Perform  this 
loop  until  the  required  number  of  measurements  have  been 
Integrated. 


Lines  3530 

Line  3550 

Lines  3570 

Lines  3600- 

Line  3630 

Lines  3650- 

Line  3670 

Line  3680 

Line  3690 

Line  3700 

Lines  3720- 

Lines  3740- 

Lines  3760- 

Line  3780 


■3540  Collect  averaged  voltage  from  voltmeter,  turn 
averaging  mode  off  and  empty  the  internal  store. 

Return  to  subroutine  two. 

-3630  SUBROUTINE  FOUR. 

•3620  Bit  7  of  the  byte  stored  in  location  59427  will  become  1 
when  a  SRQ  has  been  received. 

Return  to  subroutine  three. 

3700  SUBROUTINE  FIVE. 

K  determines  the  number  of  steps  to  be  performed  by 
stepping  motor  1. 

Z  determines  the  direction  of  steps.  The  change  in  status 
byte  acknowledges  the  receipt  of  the  instruction. 

A  null  instruction  is  sent  which  is  only  acknowledged  when 
the  previous  instruction  has  been  performed. 

Return  to  subroutine  one. 

3780  SUBROUTINE  SIX. 

-3750  Construct  instruction  containing  the  number  of  steps 

required  to  move  the  side-arm  short  circuit  to  the  right 
position  for  the  next  state. 

3770  Use  the  status  byte  and  a  null  instruction  to  determine 
successful  receipt  and  completion  of  the  instruction. 

Return  to  subroutine  one. 
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Lines  3800-3840 

Line  3820 

Line  3830 

Line  3840 

Lines  3860-4040 

Lines  3880-3890 

Line  3900 

Line  3910 

Lines  3920-3950 

Line  3960 

Lines  3970-4030 

Line  4040 


SUBROUTINE  SEVEN. 

Keep  adding  2«  radians  until  angle  is  greater  than  zero. 
Keep  subtracting  2i  radians  until  angle  is  less  than  2«. 
Return  to  main  program. 

SUBROUTINE  EIGHT. 

Construct  and  display  file  name  that  the  calibration  data 
will  be  stored  under. 

Open  channel  to  disk  drive. 

Check  if  this  filename  has  already  been  used.  If  not, 
go  to  Line  3960. 

Facility  for  overwriting  an  existing  file  or  returning  to 
main  program. 

Check  for  disk  errors. 

Store  calibration  information  and  data  on  disk. 

Return  to  main  program. 
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LISTING  OF  COMPUTER  PROGRAM 


Channel  21 

RF 

br i dge 

FORWARD  power 

meter- 

Channe 1  22 

RF 

br i dge 

REVERSE  power 

meter 

100  rem"  MULT  I  STATE  Version  2  18/1/85 

110  rent 

120  rem"  A  program  to  control  a.  MULTISTATE  REFLEC  TOMETEF:  and  solve  the 
130  rem"  non— linear  equations  involved  in  calibration  end  measur -anient 
140  rem 

150  rem"  - 

160  rem"  Channel  21  RF  bridge  FORWARD  power  meter- 
170  rem"  Channel  22  RF  bridge  REVERSE  rower  meter 

1  SO  rem  "  - - 

180  rem 

200  or en4 , 4 : openS .4.1 : open9 ,4,2:  rem"  COMMODORE  4022  pr inter 

210  open  6,6:  rem"  DATRGN  1071  voltmeter 

220  open  7,7:  rem"  RACAL-DANA  1200  switch  control  let- 

230  open  129,30s  rem"  BENTHAM  SMD-3E  stepping  motor  c)ri"e 

240  pr int#6 , "r 4t2gl " 

250  d  i  m  dt  20  >  ,  x  <  20  >  ,  y  <  20  >  ,  yr  <  20 , 5  >  ,  ro  <  20 , 5  > 

260  n2=4 j  na= 1 :  nb=4 :  z= 1 

270  input  "date  "  sda* :  input "treguencv  in  ghz.  "  :tr 
280  input  "waveguide  size  x ,  i  or  l< " ; wb* 


29© 

it 

wb*="x" 

then 

a=22 . 860 : rem " 

8.2  -  12.4 

GHz 

300 

it 

wb*  * " j " 

then 

a=15. 799 srem" 

12.4  -  18 

GHz 

310 

it 

wb*="k" 

then 

a=10.668:rem" 

18  -  26.5 

GHz 

320  it  a=0  then  print  sprint "  cannot  calibrate  that  band" : end 

330  1;  l=4tpitsgr  <  <  trl  2/299.  71 2>-<  l/<4*at2  >  > 

340  print#4 ,chr *< 1 > s "  "+wb$+  "—band  ret  lactometer " 

350  pr i nt#4  s  pr i n t#4 

360  print#4  .  "  treguencv  *  "strs"  ghz" ;spc<30> sda* spr int#4 spr int#4 
370  input  "an'-'  comments?"  sc* 

380  pr int#4 ,c*  spr int#4  sprint#4 

390  input "new  or  stored  calibration” sc*:it  cf="ne«"  then  570 
400  rem"  , - - — - — — — i 


410  rem 
420  rem 
430  til 
440  ope 
450  it 
460  i  t 
470  inp 
480  tor 
490  c  lo 
500  pri 
510  pri 
520  c*« 
530  pri 
540  tor 
550  nex 


rem"  USING  A  PREVIOUS  CALIBRATION  STORED  ON  DISK 

rem"  - - 

ti  le*="l  s "♦wb*+str*<tr  >+"ghz , sect. read" 

open2 . 8 .5  .  t  i  le*  :  rem"  COMMODORE  8050  dual  disl;  dri'-'e 

it  ds*62  then  prinf'no  calibration  data  tor  this  tr eguencv" :c lose2 send 

it  ds>28  then  print  ds* sc  lose2 send 

i  nput #2 ,  dc  t ,  me  <  1  >  ,  ms  <  2  >  .  ms  <  3  "■ ,  ms  <  4  > 

tor  i *  1  to4  s  i nput#2 . r <  i  , p <  i  >  . s <  i  >  , t <  i  ,  l<i  >  s next i 

c  lose2 

pr int#4 ,, "using  ca  libration  constants  trom  " sdc* spr int#4 
print#4,"  micrometer  " ?chr*<l>?"  rs  ps  ss  ts  Is" :prin+#4 
c*«"  99  99.99999  9.99999  9.99999  9.99999  9.9 

pr int#9 ,c* 

tor  i*l  to  4:  pr i nt#3 ,ms< i > ,r < i >  ,p< i > ,s< i >  . t <  i  >  .  1<  i  > s  print#4 
next  is  print#4s  goto  2340 


560  rem"  - 

570  rem"  NEW  CALIBRATION 
580  rem"  - 

590  input  "micrometer  settings" sms < 1 > ,ms<2> ,ms<3> ,ms<4> 

680  print#4,"  using  micrometer  settings  "s 

610  tor  i*lto3s  print#4 ,ms< i > s " , " * s  nexti i  pr int#4 ,ms<4> s  print#4 

620  rem"  - 

630  rem"  MEASUREMENTS  USING  A  LOAD  AND  SPACERS 

640  rem"  - - 

650  tor  i*l  to  3s  input  "spacer  length  <mm>"tsps  yr< i .0>*-sp*L 1 
660  gosub  3180s  rem"  SUBROUTINE  ONE 


I 


670  nz=na:ria=nb  :nb=nz  : 

680  ne:;ti 

690  a lph**yr <2.0 > -yr <  1  ,  0  >  t  beta=>T(3f0)-yr'(l  .0  > 

7O0  tor  .1  =  1  -to  4 

710  si = yr <1 , j > :  s2=vr < 2 , j  > :  s3=yr < 3 , j > 

720  c  1  =  <  s  1  -s-2  >  /  <  s  1  -s  3  > 

730  c2=  <  c.  1  ♦  <  1  -cos  <  beta'1  ))-<l-cos(»  lpha>  > 

740  c3=s  i  n< a.  lpha> -  <  c  1  its  i  n < beta  >  > 

750  theta=atri<c2/c3>  :  c2=sl/s3 
760  l  t  c3<0  then  theta=theta+pi 
770  c+=cos < theta> 

780  c3=  < c 24 cos < thata+bata> >-ct 
790  c4=2tc3/< c2—  1  > 

800  rn*  <  -c4+s-Hr  < c4tc4-4  >  >/2 
810  1  <  .j  >  *=s  1  /  <  1  +rntr  n+itmtct  > 

820  nextj 
8  30  pr  i  n-t#4 

840  a&B«bs( &lpha/pi$180>  !st=strt(a;i) 

850  b*;>abr< beta. 'pit  180>  s b$=str  t  <  ta?i  > 

860  pr  int#4, "  pos'n  0  elects  "?at;"  degs  “:b*,s"  de;is  ls"srrint#4 

£70  c$="  9  9.99999  9.99999  9.99999  9.999999" 

880  pr int#9  ,ct 

890  -for  ,j  =  l  -to  4  :pr  int#8  . .j  .vr<  1  ,.i  >  .yr  <2  ..j  >  ,yr<3...j  >.  l/i  >  snext.i 

900  r?rirvt#4  tpr  i nt#4 

910  print"  values  o-f  Is" 

920  tori=lto  4  sprint  le-t-t strt  <  l<  i  >  >  .8>  snext  i 
930  rem" 

940  rem" 

950  ren." 

960  input  "  length ( mm >  o-f  shor"t  o-f-fse-t"  so  l 

970  print  "number  o-f  settings  o-f  sliding  short"  : input  cp 
980  av=2*pi/cp :  iv=a>-.*/kl:  pt=-k  l*o l 
990  m2*="  short  offset"+strt<cil)  +  "  mm" 

1000  cS="  length  phase  pos'n  1  pos’n  2  pos'n  3  pos'n  4" 

1010  pr  int#4  ,rn2$  :pr  int#4  :pr  int#4  :p-r  int#4  :pr  int#4  .c$  :print#4 
1020  tor  i=l  to  cp 

1 030  sp= i v* < i  —  1 > :  c$=  lettS  < strS  < sp  > .6  > 

104O  print  "  set  sliding  short  to  "?c$s"  mm" 

1050  x<i)=-sp#k  l 

1060  gosub  3180:  rem"  SUBROUTINE  ONE 
1070  nz=na:na=nb  :nb=nz  :z=— z 

1080  e**="  99.999  999.9  9.99999  9.99999  9.99999  9.99999 

1090  print#9.c* 

1100  pr i nt #8 . sp , x  < i >  * 1 80/p i ; 

1110  tor  k* 1  to  4 : pr i nt#8 , yr  < i , k > : : nextk  : pr i nt#3 
1 1 20  next i 
1130  rem" 

114©  rem" 

1 1 50  rem " 

1160  tor  ic*=l  to  4  sprint "calibration  "sic 
1170  rem" 

1180  rem" 

1190  rem" 

1208  tor  i*l  to  cp:  y< i )*yr < i . ic> i  next  it  l*Kic> 

1210  tor  .j*  1  to  5 *x*x<.j >+pt :  y»y<j> 

1228  m<  j,  1  >«cos<x>  :m<  j  ,2>«-sin<x>  :m<  j  3 )«-y*co£ < x > / 1  :m<.j  .4 >=v*s i n •  >  I 
1230  m<.j  ,5>*y/  l  tnext  j 
1240  tor  ,j*l  to  3t  jl*..i+l  :.i2*.j+2 
1250  m»<  j  >  *m  <  j  1 , 1  >-m<.j2 . 1  > 

1260  mb< j>*m< j2, 1  >-tn<  1  . 1  > 


FIND  STARTING  PARAMETERS 


SOLVE  CALIBRATION  EQUATIONS  FOR  4  STATES 


MEASUREMENTS  WITH  SHORT  CIRCUITS  ATTACHED 


j  i 
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1 270  me  t  j  >  =m  t  1  ,  1  >  -m  <  j  1  . 1  > 

1280  dotj>*mtl  j>+mt  jl  .  2  >  ♦mb  <  j  >  +  rn  <  i  2 , 2  >  true  <  i  > 

1230  next  j 

1 300  al  =  <  mtl,3>^matl  >  +m  <2.3  >  ♦mb  <  1  >  «-m  <  3 . 3  >  ♦me  <  1  >  > .  ‘d>  •  <  1  > 
1310  a2= tmt  1  ,3>^ma<  1  >+m<2 ,4 >#mb<  1  >+m<3 . 4 > ♦me  <  1  >  > <  1  > 
1320  *3*<m<  1  ,3>#(ii8<  1  >+m<2,5>^mb<  1  >+m<3,5)^me<  1  )  >/dut  1  > 

1 330  a4=  <m <  1 , 3 2>+m<3,3> ♦mb <2>-t-m<4,3> ♦me < 2 > > /dw < 2 > 

1 340  a5=  <m< 1 , 4  > ♦mat' 2 >  +m <3,4  > ♦mb  < 2 >  «-m '.4.4  > ♦me < 2 >  > /dv < 2 > 
1350  a6=tm<  1  ,  5  > ♦mat 2  > +m  1 3 , 5 > ♦mb < 2 > am 1 4 , 5  > ♦me t 2 }  > /d*» C 2 > 
1360  »7*<m< 1 ,3>^ma<3>+m<4,3>^mb<3>em<5,3>^mct3>>,,d*'t3> 

1 370  aS=  t m <  1 , 4  > ♦mat  3 >  +m 1 4  , 4  > ♦mb 1 3 >  +m <5 .4) ♦me  1 3 >  > /ok-* 1 3  > 

1 380  a9=  <  m  1 1 ,  5  >  ♦mat  3  >  +m  <4,5  >  ♦mb  1 3  >  ♦m  <  5 , 5  >  ♦rue  1 3  >  >  /& 1 3  > 
139G  tl=a.3-a6:  t2=al-a4:  t3=a2-a5 
1400  t4=al-a7:  t5=a3-a.9:  t6=a2-a8 
1410  £5=  <  1 1  ♦*4— t54 12  >  / 1  1 3t:  1 4- 164:  t2  > 

1420  *4«tl/t2-*5^t3/t2 
1430  £3»a3-s4^al-s5^a2 

1440  £2=mt 1 ,5>— mt2 ,5>— s3^tm< 1 ,2>-m<2 ,2> >— s4^<m< 1 ,3>-mt 2 ,3 
1 450  £2=  t  £2— £5^ t  m <1,4 > -m  1 2 , 4  > > > / <  m 1 1 , 1 >  — m <2,1 > > 

1460  £ 1 *m 1 1 ,5>— £24mt 1,1 )— s3^m< 1 ,2>— £4tmt 1 ,3>-#5^m< 1 ,4  > 
1470  p=atn<£3/s2>  s  Hr  £2>0  -then  p=p+pi 
1480  tiQ-r-’i  giosub  3806:  p«n© 

1480  t*atnts5/*4.'>  :it  £4>0  'then  t*t+pi 
1500  n0=t:  o»o£ub  3800:  t*n0 
1518  kl*< *4 1 2+*5 1 2 > . 1  < S2T2+S3T2 > 

1 520  r-sar  <  abs  1 1  s  1  - 1  >  / 1 1  -s  1  ♦k  1  >  >  > 

1 530  £=s<nr  <  <  1  +rT2— s  1  >  /£  1  > 

1540  print  r,p,£,t:  print 
1550  ram” 

156-0  ram" 

1570  ram" 

1580  n3*ep  :n4»n2-l  :c0=5e-5  tram  convergence 
1590  tor  kc»l  to  9 

1 600  p 1 =0 : gl *0 :  ot2=0 :  r  1  *0 :  r2=0 :  r  3*=0 :  £  1  *0  :  £2=0 :  £3*0  :  s4=0 
1610  tl=0:t2»0:t3*0:t4=0: 

1620  tor  k=l  to  n3:  xb=xtk>+pt 
1630  cu*co£tp+xb> :c l=costt+xb> 

1640  tu=  1  +rtr-2^rtcu :  tl«l+s^.-2^*cl!  tm*2t\/tl 
1 650  yr«  lttu/t  1 :  dm=y t  k  >  —  vr 

1668  dm«y<k>— yr  :dl*tmttr— eu>  :d2=tm^r^s i ri t p+xb "> 

1670  d3=-2^vr/t  Wtt-e  t>  ;d4«-2^yrt*. l^£intt+xb> 

1680  pi **p  1  +dl ♦dl 
1690  ®il««tl+dl^d2 :  ct2=*i2+d2^d2 
1700  rl»rl+dl^d3:  r2*r2+d2^d3 :  r3*r3+d3^d3 
1710  *l»'sl+dl^d4 :  s2“£2+d2^d4  r  s3*£3+d3^d4  :  £4=£4+d4td4 
1720  1 1  *t  1  +d  1  ♦dm  :  t2=t2+d2^dm :  t3“t'3+d3^dm :  t4=t4+d4tdm 
1730  next  k 
1740  ram" 

1750  ram" 

1760  ram" 

1770  tl-r3^*4-*34£3s  t2»r24£4-®3^£2 :  t3«r2^*3-r3^£2 
1780  t4-rl^£4-*34£l  t  t5«r  l^*3-r3^sl  i  t6-r  1  ♦£.2-r2^£  1 
1790  ul»r2^£4-»2^*3 :  u2»«i2^r4-*2tr.2 1  u3=«i2t£3-r2^£2 
1800  u4»«tl^£4-*2^£l  t  u5*oil  ♦*3-r2^*  1  :  u6*«il^£2-«t2^£l 
1818  vl«r2^£3-*2^r3:  v2««t2^£3-*2^r2 :  o3-oi2^r3-r24r  2 
1820  y4*otl^£3-*24rl  t  o5««il^r3-r2#rl  »  v6*otl^r2-«|2^rl 
1830  t*-«i24tl-r24t2i-£2^t3i  tb-<Hl4tl-r24t4+s2^t5 
1840  tc»o»l  ♦t2-o|2^t4+£2tt6  :  td«<nl^t3-o|24t5+r24t6 
1850  ta««il  ♦t  1  -r  1  ♦  t2+*  1  ♦  t-3 :  tt-p  1  ♦  1 1  -r  1  ♦t4^-»  1  ♦tS 
1860  ta-p  1  ♦t2-otl ♦t4^-»  1  tt6 1  tb»pl^t3-c|14t5«-rl^t6 


SOLVE  4  S I MUL T ftNEOUS  EQUATIONS 


4— PARAMETER  LE0ST  SOUftRES  FITTING 
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r 


1 87© 
1 880 
1 880 
1 900 
1910 
1 920 
1 950 
1940 
1950 
1 960 
1970 
1980 
1 990 
2000 
20 1 0 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2 1 00 
2110 
2180 
2130 
214© 
2150 
2160 
2170 
2 1 80 
2190 
2200 
2210 
2220 

2248 

2250 

2260 

2270 

2280 

2290 

2300 

231© 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

243© 

2440 

2450 

2460 


d2=P  1  4te— a  1 4tb+r  1  4tc-s  1 4td 

*1-  <  +  1 4t a-+24te+434 < <t  1  *u  1  -r  1  4u2+s  1  4u3 >-+44  < 1  4v  1  -r  1  4-  -2+s  1 4-  •  3  •  • .  d. 
s2=  <  -4  1  4tb++24t+— 434 <pl  4u  1  — r  1  4u4+s  1  4u5 > ++44  < p  1  4v  1  — r  1 4v4+s  1  *5  •  • .  dl 

s3=  <41 4tc— 424toi+434<pl  4u2-«tl4u4+sl  4u6>— +44 ■>:•!  4v2— a  14'  '4+=- 1  4*  *6  1 
*4* <-41  *td+42#+h-434 <  r  1  4u3-«|l  4u5+rl4u6  >++44  <pl  4v3-al  4'  .>5+r  1  *6  >  *.  "d8 

r=r+£l  :  :  s=s+s-3:  4=++s4 

r<i=0 

•for  i  =  l  +o  n3:  a«i«x<  i  >4-p4 

ur  =  1  +r4r — 24r4cos  <  p+aa  >  :  lr=  1  +s4f-2ifs*co£ <  ++*<*> 
d+  <  i  >*v<  i  >-  I4ur/  Ir  :  r<n=ro»+d+<  i  >  t2  :  next  i 
rot=sov  < ro)/ < ri3-r>2  >  > 

i+  abs<r'<i><c0  or  abs < rfi-rm > CcO  then  2070 
pr i ri+  : for  i  =  1  ton2  :  pr i nt  s  tr-4 ■" p-3+  <  i  >  >  ; : next  i  : pr  i nt 
print"r.  m.  s.  de*  > l at  i on  =  "  :str*<roi> 
rO=rat  next  kc 
rem  " 
rern  " 
rern" 


END  OF  I TERAT I 0N8 


print  "  warn  i  ng  MM" 

print"  convergence  criterion  not  setis+i ed  " 

i+  r  -C0  then  r  =— r  :  p=p>+p  i 

nO=p :  gosuh  3800 :  p=nO 

i+  s<0  then  s=-s :  4=4+pi 

n0=+ :  oiosub  3800 :  +=n0 

d  <  1  >«o|2*  <  r34-s4-s3*s3  >  —  r  24  <  r  24s4-s-24*3  >  +s24  <  r2*s3-s2*r  3  > 
d  < 2  > -p 1 *  <  r 34s4-s34s3 > -r 1 4 • r 1 4s4-s .  1 4s3 > +s 1 *  < r 1 *f3-s 1 4r  3  > 
d  <  3  =p  1  ♦  <  oi24s4-s24s2  >  -al  4  <  «tl  4s4-s  1  4s2  >  +£1 4  <  «il  4s2-a  1 4s2  > 
d <  4  > *r>  1 4  <  ou-4r  3-r-24r  2  > -<*  1 4  < a  1 4r  3-r  1 4r  2  >  +r  1 4  < a  1 4r  2-r  1 4a2  > 
tor  il  =  l  -to  4  sfo<  i  1  .  ic>*s«tr  <abs<d<  i  1  >/d2>  >4r«i  t  next  il 
4 or  i  =  l  to  n3s  r©< i . ic>=d+<  i  > :  next  i  :  rl<ic>*r«i 
r(ic  >  «r  :  p  <  i  c  >  ssp-  :  s  <  i  c  > =s  s  +  •  i  c>  «*+ 
next  ic 

pr i n+#4 s  pr i nt#4 

pr int#4 .chr * < 1 > : "  rs  ps  ss  +s  rrns  Is" 

prin+#9  ,. "  99.99999  9.99999  9.99999  9.99999  .999999  9.99999" 

prin+#4 

+or  i*l+o4  spriri+#8  .r< i  >  .p<  i  >  .*<  i  >  .  +  <  i  >  ,rl  <  i  >  .  t<  i  > 
print #8 , b  < 1 . i > .to <2 . i > .b<  3 . i > .b<4 . i > :prin+#4 tnext i 
P-r  in+#4  .  "  residuals"  spr  in+#4 

c*= "  phase  “  :+or  i  =  1  +-c>  4:c$=c$+"  pos  'ri"+s+r  *<  i  >  inex+i 
pr-  i  n+#4  ,  c4  : pr  i  n+#4 

prin+#9."  9.9999  sr. 99999  sz. 99999  sz. 99999  sz. 99999 

tor  i  =  1  to  cp : pr i n+48  .  x < i >  : 

tor  ,j*l  to  4  :prin+#8  .r0<  i  .  i  >  :  :next.i  :pr  in+48  :rte;+i  :pr  in+#4 
input  “save  calibration  data  on  disk"  ;ct  s  it  c*0"yes"  then  2340 
aosub  3860:  ren."  SUBROUTINE  EIGHT 
rern  " 
rern" 
rerr." 

n2*2scp*4 


END  OF  CfiL  I  BRfiT  I  ON  .  STBF.’T  OF  MERSUREMENT 


pr  i  nt "  type  arid  ser.no  o+  unknown"  :  input  m2  4 
print#4 ,chrl< 13> ;chr*< 1 > » "  measurement  o+ 
•gosub  3180:  rern"  SUBROUTINE  ONE 
rem" 


STRRTING  PARAMETERS 


rem" 
rem" 

vl«v(lVKl>i  y2-v<2>/l<2>  i 
Pi-r<l>t2-yl*s<l>t2 
p2-r <  2  >  T2-v24s  <  2  >  t2 
p3-r < 3 > t2-y3*s < 3 > 12 


v3»y<3>/K3>:  41 


:m2<;  :print#4 


-vl-l»  42-V-2-1  :  43«v3-l 
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2470  «il=r<  1  >lcos<p<  1  >  >— yl#s<  1  >*cos<t<  1  >  > 

24S0  ««2=r  < 2  > ♦cos' < p < 2 >  > -v2*s  < 2 > ♦cos  < t  <2 >  > 

2490  <n3-r < 3 > ♦cos < p < 3 > > -y 34s <3> ♦cos <t<3>> 

2500  rl=yl^s<  1  > ♦s i  n  <  t  <  1  >  >-r<  1  > ♦* i n  <  p c.  i  >  > 

2510  r2“y24r*<2}^sin<t<2> >— r<2>4sin<p<2> > 

2520  r3«y3#s  <  3  >  ♦s  i  n  <  t  <  3  >  >  -r  <  3  >  ♦*  i  n<  p<  3  >  > 

2530  rem" 

2540  rem" 

2550  rent" 

2560  tl  =0i2tr3-r2tot3  :  t2*p2#r3— r2*p3 :  t3=p24i=i3-<>i24p3 
2570  dv=p  1  ♦t  1  —Hi  ♦+2+r  1  ♦  1 3 

2580  *1«<  t l*tl-t 2t Coil *r3-rl^o|3>' +t3* <<*1  ^-rl *g2> Vdv 
2590  s 2*  <  — f  1  ♦t2+t2t  < p  1  ♦r  3-r  1  ♦pS >  -t3t <p  1  *r2-r  1  4p2  >  > /dv 
2800  s3=  <  tl  *t3—f2* < p  1  ♦«t3— Hi  ♦pS > +t3t < p  1  ♦ol2-otl  *.p2 >  > /dv 
2610  um*s-’V  (  *tos  <  s  1  >  >  :  up=«tn  <  s3 /s2  >  s  i  t  s2<0  -then  up=up+p  i 
2820  n0=up :  gosub  3800  s  up=riO :  print  urn, up 
2830  rem" 

2840  rem" 

2850  retn" 

2880  n3=4 sn4=l sc0=5e-5  :r*m  convergence 
2870  tor  kc=l  to  9 
2880  pl=0  :oil=0  ;oi2=0  :t  1=0  st2=0 

2890  tor  k=l  to  n3  :  xe.=p<k>  sxb*t  <k>  srk*r<k>^um  ssk**<k>#up 
270O  cu=cos  <  ur+xe.)  s  c  l=cos  < up+xb > 

2710  tu=l+rk4rk+2^rk4cu  st l=l+sk4sk+2tsktc 1 
2720  vr= K  k  > ♦  tu/t 1 :  dm=y  < k  > — vr 

2730  dl*2/t !♦< l<k>ifcr<k>ils<rktcu>-yr**<k>ilK*k+c l>  > 

2740  d2=-2tum/t  !♦<  1  < k  >tr  < k > ♦s in< up+x*> -yrts <  1; > *s i n < up+xb >  > 

2750  p  1  =p  1  +d  1  ♦d  1  s  oil  *<*1  +d  1  *d2  s  <*2=<t2+d24d2  s  t 1 =t 1 +d 1  ♦dm  s  t2=t2+d2*dm  s 
2760  next  k 
2770  rem" 

2780  rem" 

2790  rem" 

2800  p2=<tl  sd'.»*pl^oi2“«ll^p2 

2810  um=um+  <tl ♦•h2- t2tg  1  > /dv :  up=up+  < t2tp  1  -t  1  tp2 > .‘do 
2820  roj=0stor  i  =  l  to  n3  :  pu=umtr < i > s  pl*um*s<i) 

2830  ur=  <  1  +pu#pu+2#pu*cos  <  up+p  <  i  >  >  >  /  <  1  +p  Mfp  l+2tp  Vicos  <  up+t  <  i  >  >  > 

2840  dt  <  i  )  =y  <  i  >  -  l  <  i  >  ♦ur /  Ir :  rg=r«t+dt  <  i  >  T2  s  next  i 
2850  r  «t**«ir  <  roj/  <  n3-n2  >  > 

2860  it  abs(r«t><c0  or  *tas<r0-roi><c0  then  2920 
2870  pr i nt  s  pr i nt  um , up :  pr i nt 
2880  pr  i  nt 11  r .  m .  s .  dev  i  &t  i  on  ■  *'  s  str  ♦  <  roi> 

2890  r0=roi;  next  kc 

2900  print  “  umrning  !!!!" 

2910  print"  convergence  criterion  not  setistied  " 

2920  it  um<8  then  um=-ums  up=up+pi 
293©  n0=up :  nosub  3808 :  up*n0 
2940  rem" 

2950  rem" 

2980  rem" 

2970  b  1  =*oir  <  *bs  < oi2^dv  >  >  ♦rot :  b2*so|r  <  *bs  < p  1  .'d'V  >  > ♦r ci 
2988  c»= "  Measured  estimate  res i duel  " 

2990  me-3  s.d  phese  s.d" 

3000  print#4,c# 

3010  tor  i ■ 1  to  4  sit  i>l  then  3040 

3020  c4="  99.99999  99.99999  *.99999  .999999  .999999“ 

3030  print*9.c*s  pr int#8,v< i > ,y< i >-dt< i > ,dt< i > .bl ,b2 s  goto  3660 
3040  c*-"  99.99999  99.99999  *.99999" 

3050  print#9,c*s  print#8,y< i > ,y<i >-dt< i > ,dt< i > 

3080  next! sprint#4 


MEASUREMENT  PRINTOUT 


SOLVE  2  SIMULTANEOUS  EQUATIONS 


2-PARAMETER  LEAST  SQUARES  FITTING 


SOLVE  3  SIMULTANEOUS  EQUATIONS 
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3070 

SOSO 

3030 

310O 

3110 

3120 

3130 

3140 

3150 

3180 

3170 

3180 

3130 

320O 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3230 

3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3330 

3400 

341© 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3430 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

357© 

3580 

3530 

3600 

3610 

3620 

3630 

3640 

3650 

3660 


r  $=str t  <  urn  >  :  r*=s+r *  •;  up*  1  80.-  ‘p  i  > 

print#4  ;pr  irvt#4 ,  "rho  ot  ":rn2*;"  =  "srts"  at  "  zr-  t  s  "  decjr  eef" 
nz*nasna=nb  :nb=nz  sz=-z 

i nput " imo+her  deuice"  sc*  : it  c*="yes"  ■then  2340 
c  lc<se4  :  c  lose6 :  c  lose 7 s  c  loseS :  c  loseS  s  c  lose  1 23 
end 
rew  " 
rent" 
ren," 


END  OF  PROGRAM 


rern 
rent " 
rent" 
rent " 

w2-l:«o*ub  3330 srem"  SUBROUTINE  TWO 
z  1  *•...•  1  :  z2=u2  :  «>2=0 

z*="  i  "  :  ojosub  3650  srem"  SUBROUTINE  FIVE 

•tor  ,i*na  "to  nb  step  z : ojosub  3330 srem"  SUBROUTINE  TWO 

y  1  <  j  >  =u  1  :  m2  <  j  >  »u2 

it  ,j*rib  then  3230 

rns=ms  <  j +z  >  -ms  <  ..i  > 

z*="d"s  it  r..sC0  then  z*="i" 

»osub  3720  s  rern"  SUBROUTINE  SIX 
next  j 

z*="d" s  ojosub  3650  srem"  SUBROUTINE  FIVE 
u<2=l  s-aosub  3330 srem"  SUBROUTINE  TWO 
z4=vl  :  z5='v'2  s  ui2=0 

ac=l sfor  j»na  to  nb  step  z:bc*l-ac 
za=  <  actz  1  > t < bc*z4  >  s  zb=  <  actz2 > ♦  <  bc*z5  > 

y  <  j  >  =  <  zb  1 2-m2  < ..i  >  t2 > / < zat  2-u  1  <  j  >  t 2  >  : y*=  lett*  < str  t  <  v  <  ..i  >  >  .7  > 
printy* sac«ac-< l/3>  syr<i  >«y<j>  snextj 
return 
rem" 
rent " 
rent" 

print#7, "rx" 
nd=8  sit  w2= 1  then  nd=4 

print#7,"c21x"  sojosub  343©  srem"  SUBROUTINE  THREE 
ul=ual<ul > 

pr  i  nt#7  ,  "  c22x  "  s  gosub  3430  s  rern •'  SUBROUT  I NE  THREE 

m2 I  <  M*  > 

return 


SUBROUTINE  TWO 


SUBROUTINE  ONE 


re  m  " 
rem  " 
rern  " 

print#6, "al " storaa=l to 100 snextaa 

tor  k«l  to  nd s pr i nt#6 ,  " ffl "  s ojosub  3570  s  rern"  SUBROUTINE  FOUR 
nextk  s  i  nput#6 ,  u*  ,  ui* 
pr  i  r.t#6  ,  "  *0 10 "  spr  i  nt#7  ,  "rx " 
return 
rern" 
rem" 
rem" 

a«peek< 59426 > 
tor  zx*l  to  5000 

b*i»eek < 53427 >  and  128 sit  b-0  then  next  zx 
it  zx>4938  then  print"sr«i  not  detected" 
return 
rem" 
rem" 
rem" 


SUBROUTINE  FIVE 


SUBROUTINE  FOUR 


SUBROUTINE  THREE 
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3670  k*="0200w" 

3680  pr  i  nt#  1 29 r*+k*  sit  stO0  -then  3680 
3690  pr  i  nt#  1 29  ,  "  i  0OO0v  “  s  i  t  rtO0  then  3690 
3700  return 
3710  rem" 

3720  rem  “ 

3730  rem " 

37 40  k/i*420tms  s  k*  len < str#  <  k*i > >  :  kl«m i  d*  <  str*  K  kX  >  ,  2 ,  k— 
3750  it  len<k*><4  then  k**"0"+k*s  «oto  3758 
3760  print#129,z*+k*-t-"v"  sit  stO0  then  3760 
377©  pr  int#129,*'i 0000v"  sit  *tO0  then  3770 
3780  return 
379©  rem" 

3808  rem" 

3818  rem" 

3820  it'  n0<8  then  n0=n0+2*pi  soioto  3820 
3838  it  ri0>2#pi  then  n0=n0-2*pi  :<aoto  3830 
3840  return 
3858  rem " 

3860  rem" 

3870  rem" 

3888  cr*=chr*<:i3>  sti  le*="l  s “+wb*+str*< tr>+"3hz" 

3890  print"  +  i  le  name  will  be  “ sti le* sprint 
3908  open2,8,5  ,ti  le*+"  ,se-=i, write" 

3910  it  drO€3  then  3960 

3920  print"  a  ti le  ot  this  name  already  exists" 

3930  input "do  you  wish  to  over-write  it" sc* 

3940  it  c*="no"  then  pr int#2, sc  lose2 * return 
3956  print#2,">" sti le*="G"+ti le* sc lose2 :«oto  3900 
3960  it  ds>20  then  printds* sc lose2 s return 
3970  ml=ms<l>:  m2*ms<2>  s  m3*ms<3>:  m4*ms<4> 

3980  pr i nt#2 , da* ; cr *  s  m 1 * cr*  s  m2 s  cr * s m3 ; cr *  t m4 scr* s 
3990  tori=lto4 

4000  dl»r< i >  sd2=p< i >  sd3*s< i >  sd4=t< i >  sd5=  Ki) 

4010  print#2,dl scr* sd2 scr* ;d3 scr* sd4 scr* sdSscr*  s 
4020  nextisclose2 

4030  print  sprint"  tile  "stile*s"  written  to  disk" 

4040  return 
ready . 


SUBROUTINE  EIGHT 


SUBROUTINE  SEVEN 


SUBROUTINE  SIX 
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CONCLUSION 


The  equipment  described  is  very  versatile,  easy  to  use  and  easily  adapted 
to  any  waveguide  band.  It  produces  excellent  results  at  a  relatively  low  cost. 
The  accuracy  is  dependent  upon  the  accuracy  of  the  standards  used  to  calibrate 
the  ref lectometer  and  the  repeatability  of  the  side-arm  short  circuit.  The 
repeatability  can  be  very  tightly  controlled  if  a  rotating  short  circuit  is 
used.  The  use  of  a  non-contacting  short  circuit  degrades  the  performance 
slightly.  Other,  faster,  methods  of  solving  the  equations  can  easily  be 
applied  and  if  an  electronic  method  of  changing  the  state  repeatability  could 
be  found  this  would  radically  improve  the  speed  and  usefulness  of  the  equipment. 
The  principle  can  be  applied  to  media  other  than  waveguide,  eg  coaxial  or 
stripline,  where  suitable  standards  can  be  found. 
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Line  2660 


Set  number  of  data  point*  and  criterion  for  aucceasful  convergence. 


Line  2670  Perform  up  to  nine  Iteration*. 
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